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Abstract. The HEGRA system of four Imaging Atmo- 
spheric Cherenkov Telescopes (IACTs) has been used to 
extensively observe extragalactic objects. In this paper 
we describe the search for TeV emission from nine very 
promising potential TeV sources, namely eight "high fre- 
quency" BL Lac objects (HBLs), and the object "BL Lac- 
ertae" itself. These objects were observed during 1997 and 
1998 seasons, with total integration times ranging between 
one and fifteen hours. No evidence for emission was found 
from any of these objects and the upper limits on the in- 
tegral energy flux above ~750 GeV are on the level of 
a few times 10 -12 erg cm~ 2 s~ 4 . For the two objects BL 
Lacertae and 1ES 2344+51.4, we discuss the astrophysical 
implications of the TeV flux upper limit, using also infor- 
mation from the X-ray and 7-ray bands as measured with 
the All Sky Monitor (ASM) of RXTE (1.3-12.0 keV) and 
with EGRET (30 MeV - 20 GeV). 



Send offprint requests to: C. Renault (rcecile@in2p3.fr) 
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1. Introduction 

The phenomenology of BL Lac objects has led to a clas- 
sification of these sources into two sub-populations: high- 
frequency BL Lac objects (HBLs) and low-frequency ones 
(LBLs) (see e.g. Fossati et al. 1998). The HBL objects have 
synchrotron and Compton peaks at relatively high ener- 
gies (Ghisellini et al. 1998), and thus are good candidates 
for TeV emission. With the HEGRA system of IACTs the 
significant detection (e.g. > 3 a) of a flux comparable to 
the flux from the Crab nebula takes less than 15 minutes 
for sources close to zenith position. The BL Lac observa- 
tions presented in this paper are thus sufficient to probe 
energy fluxes down to the level of 10~ 12 erg cm -2 s~ 4 
(> 750 GeV). 

The two BL Lac objects Mkn 501 and Mkn 421 are 
well established and well studied TeV sources. Mkn 501 
showed during 1997 spectacular outbursts with diurnal 
flux levels reaching 10 Crab units (e.g., Aharonian et al. 
1999a, Samuelson et al. 1998, Djannati et al. 1999) and 
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with a spectrum extending up to at least 16 TeV (Aharo- 
nian et al. 1999b). For a better understanding of the intrin- 
sic properties of this class of sources and for determining 
the amount of intergalactic extinction of the TeV photons 
due to pair production processes with the Cosmic Infrared 
Background radiation (CIB) photons, the detection of fur- 
ther BL Lac objects is of utmost importance. Thus the 
search and the detailed study of new extragalactic TeV 
sources is one of the major objectives of the HEGRA ex- 
periment. 

The paper is organized as follows: We describe the data 
samples used in the analysis in Section 2 and the method 
used for computing flux upper limits in Section 3. The ex- 
tinction by the CIB is discussed in Section 4, and in Sec- 
tion 5 we give the experimental results. In Section 6 and 
7 we focus on two objects: (1) 1ES 2344+51.4 detected 
during December 1995 by the WHIPPLE Chcrcnkov tele- 
scope (Catanesc et al. 1997), and BL Lacertae which, in 
July 1997, showed a strong flare detected in gamma-ray 
by EGRET, in the visible (Bloom et al. 1997; Madejski et 
al. 1999) and in X-ray wavelengths. A strategy to optimize 
the search for new extra-galactic sources is briefly outlined 
in Section 8, and the conclusions are given in Section 9. 

2. Data 

The HEGRA experiment, located on the Canary Island 
La Palma at the Observatorio del Roque de los Mucha- 
chos (2200 m a.s.L, 28.75°N, 17.89°W), consists of several 
arrays of particle and Cherenkov-detectors dedicated to 
cosmic ray research (Lindner 1997, Barrio et al. 1998). 

Before September 1998, the CT-System consisted of 
four, (now it consists of five) telescopes with 8.5 m 2 mir- 
ror area each. Each telescope is equipped with a 271 pixel 
camera with a pixel size of 0.25°, covering a field of view 
of 4.3°. The cameras are read out by 8 bit 120 Mhz 
Flash-ADC systems. Details about the CT-system and 
the performance of the stereoscopic air shower observation 
method are given by (Daum et al. 1997) and (Aharonian 
et al. 1999a). 

Table 1 summarizes the observations and redshifts of 
the nine BL Lac objects considered in this work. The first 
7 objects with right ascension < 9 h have been observed 
during October and November 1997. The data have been 
taken with only 3 telescopes in the system because the 
telescope CT4 was not operational from October 16 to 
November 15 1997. The last two sources have been ob- 
served with the 4-CT system, BL Lacertae in July/August 
1997 and May/ June 1998 and 1ES 2344+51.4 in Decem- 
ber 1997. 

3. Upper limits in Crab units and in flux units 

In the following analysis, measured fluxes or flux upper 
limits are first determined in Crab units and are only 
subsequently converted into absolute flux values or flux 
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Table 1. Summary of the data samples for the 9 BL Lac ob- 
jects. (9) denotes the mean zenith angle of the observations 
and n.t. the number of telescopes in the stereoscopic system. 



upper limits. Results in Crab units have the advantage 
of relying exclusively on measured data and of being free 
from systematic errors due to the Monte Carlo simula- 
tions of the air showers and the detector response. More- 
over, as several Bl Lac objects have been observed with 
only three telescopes in the CT-system, we can directly 
compare their observations to Crab data taken under the 
same conditions, at the same epoch, avoiding thus spe- 
cific simulations. The trigger rate is constant within 5% 
for a given zenith angle (Aharonian et al. 1999a), thus no 
strictly simultaneous data can be compared safely. 

For each source and each zenith angle interval I za the 
number of events in the ON-source region (ON s ), the num- 
ber of events in the OFF-source region (OFF s ), and the 
observation time T s is determined. To maximize the statis- 
tics the analysis is based on "loose" cuts: the mean scaled 
width of the showers (Konopclko 1995; Daum et al. 1997) 
has to be smaller than 1.2 (to retain w 80% of photon in- 
duced showers) and the squared angular distance of the 
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reconstructed shower direction from the source direction 
has to be smaller than 0.05 deg 2 . For each time period 
of fixed experimental conditions, a reference Crab data 
sample was analyzed (compare Table 2) and the numbers 
ON c , OFF c , and T c were determined for all zenith angle 
intervals. As Crab observations with three telescopes have 
been performed only for several days, the flux upper lim- 
its derived with this reference data are considerably higher 
than for the data with 4 telescopes in the system. 

Using the probability density function of the number of 
source events, we compute the upper limit of the number 
of counts MAX S from T s hours of source observations at 
99% confidence level (Helene 1983). Similarly, we calculate 
the lower limit MIN C of the number of counts for T c hours 
of Crab observations. We compute the upper limit in Crab 
units UL S from: UL S = x ^. 

The energy threshold E t h is computed by Monte-Carlo 
for each source as a function of the mean zenith angle of 
the object during the observations. For the CT-system, 
the energy threshold scales with zenith angle 8 roughly as 
cos~ 25 (#) (Konopelko et al. 1998). Only data with good 
weather conditions are used and phototube voltage fluc- 
tuations are corrected. 



MJD of the Crab observations 


hours 


3 telescopes 


50747 to 50750, 50760, 50777, 50778 


15h40 


4 telescopes 


50721, 50731, 50732, 50788, 50789, 50806, 50831 




50809 to 50811, 50835 to 50837, 50864 to 50866 


54h45 


50869, 50872 to 50874, 50891 to 50893, 50902 





Table 2. Information about the observations of the Crab Neb- 
ula used for computing flux upper limits in Crab units. 



Assuming a source energy spectrum, the conversion of 
upper limits in Crab units into upper limits in absolute 
flux units is straightforward. In the following we use two 
slopes for observed source spectrum: dN/dE oc E~ 26 , as 
measured for the Crab nebulae around 1 TeV (Konopelko 
et al. 1998) and a steeper one dN/dE oc E~ 3S . Above 0.5 
TeV, the integral Crab flux is F-2.6— 5 ■ 10~ n cm _2 s _1 . 
In the other case, we normalise the flux in order to get the 
same integral flux in this range. 

The upper limits on the integral flux are then computed 
from: 

= UL S x 1.7 ■ 10- 11 E~^ 6 cm" 2 s" 1 , and 
UL_ 3 , e = UL S x 8.5 ■ 10~ 12 E~ h 2 - 6 cm" 2 s" 1 . 

4. Correction for intergalactic absorption 

Interaction of TeV gamma-rays with intergalactic infrared 
photons by pair production processes probably modifies 
substantially the intrinsic TeV spectrum emitted by the 



source (Nikishov 1962; Gould & Schreder 1965; Stecker et 
al. 1992). The most important domain for us is 1-10 /.tm as 
these photons can interact with TeV gamma-rays around 
the HEGRA threshold. The curve in Fig. 1 presents up- 
per limits and measurements in relative agreement with 
recent modelling of CIB (Malkan & Stecker 1998; Primack 
et al. 1999). Gamma photons below 2 TeV effectively in- 
teract only with CIB photons below few /xm, so only Pri- 
mack models can be used. As the purpose of this section is 
to confirm ideas, we choose the LCDM cosmology model 
(cosmological constant + cold dark matter) which gives 
higher density and is therefore more conservative than the 
HCDM one (hot+cold dark matter). 

For a non-evolving CIB density according to this 
model, the amount of absorption for 0.1 to 20 TeV photons 
is shown in Fig. 2 for three of the studied BL Lac objects. 
Strictly speaking, the CIB energy density depends on red- 
shift due to evolutionary effects, which are not well known. 
However, the redshift dependence could be neglected for 
redshifts lower than 0.15 since we do not expect significant 
evolution on very short time scales. 

Between 500 GeV and 1 TeV, 25% of the flux is ab- 
sorbed even for a source as close as Mkn 501 with z=0.034; 
for the source 2E0829. 1+0439 with a redshift of 0.18, this 
percentage rises to 75%. The figure shows that the CIB 
absorption becomes very large for objects with redshifts 
above ~ 0.15. 




10" 1 1 10 10 2 A,(^m) 

Fig. 1. Experimental constraints on the CIB energy density to- 
gether with LCDM model by Primack et al. 1999 (full line) and 
the maximum CIB given by Malkan & Stecker 1998 (dashed 
line). The triangles show ISO measurements at mid-infra-red 
(Stanev et al. 1998); the open circle shows a tentative detection 
at 3.5/im (Dwek et al. 1998a); all other measurements, lower 
and upper limits are from a recent compilation by (Dwek et 
al., 1998b). 
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Fig. 2. Percentage of VHE photons reaching the observer after 
absorption by the CIB using LCDM model by Primack et al. 
1999. 



5. Experimental results 

The observations described in this paper do not reveal pos- 
itive evidence for VHE emission from any of the nine stud- 
ied sources. Table 3 summarizes the upper limits in Crab 
units UL S , the mean energy thresholds of the observations, 
and the upper limits on the integral fluxes, assuming the 
two different slopes for the spectrum, above the energy 
threshold with and without correction of the absorption 
by the CIB (Primack LCDM model). The threshold be- 
ing defined as the product of the acceptance of the de- 
tector and the source gamma-ray spectrum, it effectively 
depends on the slope of the gamma-ray flux. So the table 
give thresholds for the Crab slope (-2.6), and an effective 
threshold 30% lower is used for computation of the steeper 
model. 

All upper limits have a confidence level of 99%. As 
mentioned above, the flux upper limits in Crab units are 
free from systematic errors. The major systematic uncer- 
tainty of the upper limits in absolute flux units derives 
from a 15% uncertainty in the energy scale of the CT- 
System. This 15% uncertainty translates into a flux un- 
certainty of approximately 20% and 40% for a integral 
source spectral index of -1.6 and -2.6 respectively. 

Upper limits of a few times 10 -12 cm~ 2 s _1 are obtained 
after a few hours of observation. The correction of the 
absorption by the CIB of Fig. 1 (computed up to 3 TeV) 
increases the upper limits by a factor of 2 to 60 (for a 
slope of -2.6) depending on the effective threshold and the 
source redshift. 

Of course the effect of the CIB absorption makes the 
observed spectrum steeper and then the impact of a very 
hard intrinsic spectrum is attenuated (upper limits for the 
two slopes are less different after absortion than before). 
We note that the upper limits are derived from assumed 
slopes for spectra and a model for the CIB intensity. These 
uncertainties make upper limits rather qualitative. 
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Table 3. 99% confidence level upper limits derived in this 
work: (a) right ascension of the source, 2200a designates 
BL Lacertae in 1997, and 2200b BL Lacertae in 1998, (b) mean 
energy threshold of the observations in GeV for the -2.6 slope 
(the effective threshold is approximately 200 GeV lower for 
the -3.6 slope), (c) upper limits in Crab units, (d) upper lim- 
its on the integral flux above the threshold energy in units of 
10~ 12 cm~ 2 s~\ and (e) upper limits on the integral flux be- 
tween the threshold energy and 3 TeV after correction for the 
CIB absorption (using LCDM model by Primack et al 1999) in 
units of 10~ 12 cm~ 2 s _1 . The slope of the assumed differential 
spectrum is -2.6 or -3.6 (see text for details). 



6. BL Lacertae (1ES 2200+42.0) 

Several experiments extensively observed this source fol- 
lowing the detection of a strong flare by EGRET in the 
7-ray regime (only several hours observation time) and 
simultaneously in the optical regime around July 19th, 
1997 (MJD 50648) (Bloom et al. 1997; Madejski et al. 
1999). Because this flare occurred during a full moon pe- 
riod, the TeV observations started only about ten days 
after the detection of the flare. Our observations yielded a 
flux upper limit of 4.5-10" 12 ph cm" 2 s" 1 (E >580 GeV). 
The CAT group obtained a similar upper limit: 11.4 • 
10- 12 phem- 2 s- 1 («11% of the Crab) (E >300 GeV) 
(Barrau, private communication). Earlier measurements 
at VHE energies were performed in 1995 with the Whip- 
ple telescope following a weak flare detected by EGRET 
with an integral 7-ray flux of 4 ± 1.2 • 10~ 7 ph cm~ 2 s _1 
(E >100 McV), i.e. three times lower than the flux of the 
July, 1997 flare. The Whipple observations (40 h) yielded a 
flux upper limit of 5.3 10~ 12 ph cm" 2 s" 1 (E >350 GeV) 
(Catanese et al. 1997). The results from the three VHE 
experiments HEGRA, CAT, and Whipple are shown in 
Fig. 4. 

Information about the X-ray activity between Febru- 
ary 1996 and August 1998 in the energy region 2-12 keV 
is provided by the All Sky Monitor (ASM) on board the 
Rossi X-Ray Timing Explorer (RXTE) (Remillard 1997). 
We determined the ASM count rates from the "definitive" 
results obtained through analysis of the processed data by 
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the RXTE ASM team at MIT; data have a dwell duration 
larger than 30 seconds and a flux fit with a reduced x 2 ~ 
value below 1.5. The light curve is extracted using the 
"ftools 4.0" package. Unfortunately, the mean count rates 
are very low, i.e. typically around 10 _2 -10 _1 Hz and the 
observation frequency is less 1 "run" of 90 seconds per 
hour, making it difficult to get information about the 
source activity on time scales as short as one day. So data 
are binned in 5 week bins, yielding flux estimates with 
acceptable statistical errors. 

The light curve is shown in Fig. 3. A pronounced lu- 
minosity increase during June 1997 followed by a slow de- 
crease during the following 5 months can be recognized. 
HEGRA data were taken in July/August 1997 when the 
source was still active. In 1998, Bl Lacertae was in a low 
X-ray activity. 



iL Lac 



+ 



600 800 
Time (MJD-50000) 



Fig. 3. Light curves computed from data of the ASM/RXTE 
detector (1.3-12.1 keV) with a binning of five weeks. Observa- 
tions range from February 1996 to August 1998. 



Note that EGRET was able to detect the source, and 
the TeV instruments were not. Intriguingly in the case of 
Mkn 501 the opposite happened: the source was bright in 
the TeV energy regime but could hardly be detected at 
MeV/GeV energies. 

Model calculations by Bottcher & Bloom (1998) and 
Madejski (1999) predict TeV flux significantly higher than 
the upper limits obtained by HEGRA (and also by CAT at 
the same period). We explain the non-detection at VHE 
energies by two reasons. First, the absorption by a CIB 
density as in Fig. 1 reduces the observable flux by about 
50%. Secondly, the HEGRA measurements were taken 
roughly 10 days after the detection of the GeV-flare. Tak- 
ing into account the strong variability observed in the two 
BL Lac objects Mkn 501 and Mkn 421, the VHE flux of 
BL Lacertae could have easily decreased by a factor of >5 
in 10 days. 

7. 1ES 2344+51.4 

The spectral energy distribution of 1ES 2344+51.4 is 
shown in Fig. 6. This source was observed by the Whipple 
group during the winter months of 1995/96 and 1996/97 
(Catanese et al. 1998). Whipple detected a flare on Decem- 
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Fig. 4. Spectral energy distribution of the object BL Lacertae. 
Data are taken from (Perlman et al. 1996; Webb et al. 1998; 
Catanese et al. 1997); different measurements well before and 
after the 1995 and the 1997 flare (circles), measurements 
around the 1995 flare (squares) and measurements around the 
1997 flare (triangles). HEGRA CT-system data are represented 
by filled symbols while others are represented by open ones. 
The black symbols indicate the raw limit, the grey ones the 
limit after deconvolution by the CIB absorption (LCDM model 
by Primack et al. 1999). 



ber 20 th 1995 with a flux of 6.6±1.9 lO" 11 ph cm~ 2 s" 1 
(~ 63% of the Crab); the mean flux (excluding the 
flare emission) during the first winter was estimated as 
1.1±0.4 lO" 11 ph cm" 2 s^ 1 («11% of the Crab). For the 
period 96-97 an upper limit of 8.2 10~ 12 ph cm -2 s _1 
(«8% of the Crab) was reported (all integral fluxes above 
350 GeV). 

The flux limit computed from the 1997/1998 HEGRA 
CT-system data is 2.9 10" 12 ph cm" 2 s" 1 (E > 1 TeV). 
The smooth curve in the 2-10 keV band (Fig. 5) indicates 
that the source has been in a low and stable state for 
1.5 years. The source was not monitored by the ASM dur- 
ing the winter 1995/1996, therefore, unfortunately, the 2- 
10 keV activity during the state of VHE-emission could 
not be examined. 



1ES 2344+51.4 



Time (MJD 



900 
50000) 



Fig. 5. Light curves computed from data of the ASM/RXTE 
detector (1.3-12.1 keV) with a binning of five weeks. Observa- 
tions range from June 1997 to August 1998. 
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alert-system which allows us to react to a flare detection 
in any chosen wavelength within several hours. 

Fig. 7 shows the observation time required to achieve 
a detection for sources with redshifts ranging from 0. to 
0.2. The intrinsic emission is assumed to be 10% to 100% 
of -F-2.6 (assuming a slope of 2.6) or of F_ 3 . 6 (with a slope 
of 3.6, see Section 3). It can be recognized that, depending 
on the redshift of the source, observation times of several 
hours could suffice for a significant detection. 



v (Hz) 



Fig. 6. Spectral energy distribution 

of the source 1ES 2344+51.4. Data are taken from Catanese 
et al. 1998 and Qin and Xie 1997: measurements of different 
periods before and after the 1995/96 winter period (circles), 
mean flux value of the high state during the 1995/1996 win- 
ter period (square) and flux during the 1 day flare observed 
by Whipple (triangle). The HEGRA CT-system upper limit is 
shown as a filled symbol while complementary data are rep- 
resented by open symbols. The black dot indicates the raw 
limit, the grey one the limit after deconvolution by the CIB 
absorption (LCDM model by Primack et al. 1999). 



8. Observational strategy 

It is widely believed, that TeV gamma-rays of BL Lac ob- 
jects are produced by ultra-relativistic electrons, which 
emit synchrotron radiation in the keV band and produce 
the TeV photons due to Inverse Compton scattering of 
soft target photons (e.g., Ghisellini et al. 1998). The X- 
ray activity of a source gives crucial information about the 
presence of high energy electrons which could also pro- 
duce TeV radiation. A major uncertainty arises since the 
keV photons could either be synchrotron photons of very 
high energy electrons (the actual energy of the electrons 
certainly depends on the value of the magnetic field) or 
could be inverse Compton photons produced by electrons 
with moderate energies. In the latter case, TeV emission 
of the source would be less probable. 

The ASM/RXTE data provide a unique tool for moni- 
toring the light curve of many X-BL Lac objects, but with 
a moderate sensitivity. In the case of Mkn 501, the state 
of increased activity lasted several months (Aharonian et 
al. 1999a). The ASM data could be used to detect such 
a state by averaging the data over several days or weeks. 
Such a large bining time does not allow us to search for 
short-term variability but is dictated by the low sensitiv- 
ity of the detector. Moreover, as for Mkn 421, the range 
2-12 keV can be located around the "pivot point" of the 
synchrotron spectrum (the hardening starts at the energy 
of the pivot point): it is not a completely safe indicator of 
activity, but it is certainly the best one available. When- 
ever such a state was detected, observations could be ini- 
tiated. Such a strategy could be complemented with an 
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Fig. 7. Expected observation time for achieving a 3-<r detection 
for a flux (before absorption by the CIB) of 10%, 20%, 50% or 
100% of F-2.6 (full line) or of F-3.6 (dashed line). Features of 
the HEGRA 4 CT-system and an energy range of 0.5-3 TeV 
have been assumed for the calculations; the results are valid 
for observations with zenith angles below ~20°. 



9. Conclusions 

In this paper we report upper limits for a sample of eight 
X-ray selected BL Lac objects and the object "BL Lacer- 
tae" observed in 1997/1998 with the HEGRA CT-system. 
Seven sources were studied with three telescopes during 
one to ten hours. The TeV upper limits range from 15% 
to 65% of the Crab flux. The conversion in units of cm~ 2 
s _1 is done assuming the slope of the spectrum and a 
model of CIB absorption. 

The sources BL Lacertae and 1ES 2344+51.4 were ob- 
served with the 4-CT-system for approximately 15 h each. 
BL Lacertae was flaring in July, 1997 in optical and soft 
gamma-ray. The source 1ES 2344+51.4 is a very good 
VHE candidate due to its similarity to the well established 
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strong TeV sources Mkn 421 and Mkn 501. Our upper lim- 
its are at the level of ~10% of the Crab flux. 

The HEGRA IACT system could detect TeV sources 
out to a redshift of 0.15. Beyond this, intergalactic absorp- 
tion is expected to reduce the >500 GeV flux consider- 
ably. Since BL Lac sources are known to be very variable, 
a search for TeV emission is much more promising when 
guided by observations in other wavelengths, especially by 
observations in the X-ray energy band. 
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